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REAL TIME OIL RESERVOIR EVALUATION 
USING NANOTECHNOLOGY 

RELATION TO OTHER APPLICATIONS 

5 

This application is a Continuation In Part of U.S. Ser. No. 
11/178,079, filed 8 Jul. 2005, and of U.S. Ser. No. 11/416, 
505, filed 28 Apr. 2006. 

ORIGIN OF THE INVENTION io 

This invention was made by one or more employees of the 
U.S. government. The U.S. government has the right to make, 
use and/or sell the invention described herein without pay- 
ment of compensation, including but not limited to payment 15 
of royalties. 

FIELD OF THE INVENTION 

This invention relates to use of nanotechnology-based 20 
chemical sensors to evaluate chemical and physical param- 
eters associated with oil or gas well exploration, development 
and operation. 

BACKGROUND OF THE INVENTION 25 

Current mechanical, electromagnetic or optical sensor 
technologies are relatively low-tech and passive, and do not 
acquire data at sufficient distances to permit reservoir man- 
agers to fully comprehend the chemical composition, volume 30 
and dynamics of the petroleum in a given reservoir. The most 
advanced tool for monitoring the wells is a multiphase flow 
meter. It monitors the flow rate of oil, water and gas. There is 
no tool or methodology that can provide the information 
about the quality of the oil such as hydrocarbon content and 35 
the oil to water ratio, which related to the volume and pro- 
duction enhancement methods (i.e. steam injection). 

In large part, the future of the oil and gas industry depends 
on the ability to better understand the volume and dynamics 
of a reservoir to optimize production and avoid damaging the 40 
reservoir or interrupting flow through over-production or 
other production enhancement methods such as steam injec- 
tion. The disclosed chemical nanosensor network, combined 
with the physical (micro)sensors such as humidity, tempera- 
ture and pressure, acoustic or electromagnetic wave, form a 45 
monitoring system can substantially improve the quality and 
production of the oil by monitoring the chemical composition 
in wells and then feed back the information to the decision 
maker to modify and fine tune the production enhancement 
methods in real time to improve the oil quality and control the 50 
quantity. Without adequate and reasonably complete data for 
a candidate wellsite, the success rate is presently 20-25 per- 
cent and has not increased much in the last 20 years. 

The geographic market for this technology is global. Large 
and small U.S . and North American reservoirs are applicable, 55 
so as those worldwide reservoirs. The potential of this new 
technology to dramatically increase oil production will have 
a significant downward impact on world prices. Even a mar- 
ginal enhancement of the sensor capability of oil and gas 
monitoring will produce exponential benefits. The potential 60 
of heavy oil is entirely technology driven. In situ production 
methods (as opposed to open pit mining methods) require 
steam injection that causes the liquefaction of tar sands 
deposits . But without solid data regarding the dynamics of the 
geological system, production techniques result in only a 65 
20-25% recovery rate. The future ability of conventional oil 
and gas production to meet world energy demand while also 
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reducing political tensions depends almost entirely on 
increasing the productivity of known reserves through new 
technology significantly upward from the current average of 
approximately 33%. There currently is no active monitoring 
system on the market capable of meeting that requirement. 

What is needed is an integrated system and method for 
estimating one or more qualitative or quantitative parameters 
associated with an underground reservoir of a fluid mineral 
(e.g., oil or natural gas) that permits a more accurate assess- 
ment of the economic potential of the reservoir. Preferably, 
this assessment should include one or more of local (under- 
ground) values relative humidity, temperature, gas pressure, 
fluid level, and/or presence/absence of one or more target 
molecules, such as C m H„. Preferably, this assessment should 
also permit an estimate of direction and flow rate (represented 
by a current vector J) of the fluid in response to pressurization 
or other perturbation of the fluid mineral resource. Preferably, 
this approach should apply to evaluation of a “new” well 
(identified but not yet developed), of a producing well and of 
a capped well. 

SUMMARY OF THE INVENTION 

These needs are met by the invention, which provides a 
system and associated method for receiving and evaluating 
relevant parameter values, in exploration and operation of an 
oil well or a field of adjacent oil wells. One or more nano- 
structures (“NSs”), including but not limited to carbon nano- 
structures, in an NS array is coated or doped (“loaded”) with 
a selected chemical substance (e.g, chlorosulfonated polyeth- 
ylene, hydroxypropyl cellulose or mono -layer protected clus- 
ters (“MPCs”) of Au. The chemical substance is chosen so 
that, when an environment with a selected environmental 
parameter value (e.g., relative humidity, temperature, fluid 
pressure, fluid motion, presence/absence of a target mol- 
ecule) is present, a selected electrical parameter (e.g., current, 
voltage difference, resistance, conductance) will have a cor- 
responding measured parameter value, which preferably var- 
ies monotonically with the electrical parameter value. Alter- 
natively, measurement of a single electrical parameter value 
may correspond to values of each of a group of distinct 
environmental parameters. 

This technology will provide two-dimensional fluid mea- 
surements by monitoring the chemical composition in gas or 
liquid phase and the physical properties of oil level in the 
wells to have comprehensive digital information from each 
well that distributed in the reservoir with a horizontal extent 
of around 1-10 miles and 2-2000 feet in depth. Because the 
chemical and physical measurement are digital, satellite com- 
munications technology can be utilized to bridge the distance 
between people and data. Reservoir data can be centralized to 
allow an oil and/or gas company and its clients to monitor and 
analyze underbalanced drilling projects without traveling to 
the wellsite. This technology can reduce risk, cut the amount 
of time needed for carrying out exploratory operations, and 
reduce the size of a reporting system placed in a wellsite 
down-hole. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 illustrates an array of nano sensors, each using a 
network of SWCNTs to connect two electrodes according to 
an embodiment of the invention. 

FIGS. 2A and 2B illustrate an interdigitated electrode 
(IDE) array (2A) and an SEM image of a Pd-loaded section of 
the IDE according to an embodiment. 
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FIG. 3 graphically illustrates variation of a measured elec- 
trical parameter with time for different concentrations. 

FIGS. 4 A and 4B compare relative response of SWCNTs 
coated with chlorosulfonated polyethylene or with hydrox- 
ypropyl cellulose for presence of N0 2 , Cl 2 , HC1, and CH 4 . 

FIG. 5 is a flow chart for practicing the invention using one 
or more selected coatings and dopants (e.g., Pd) to detect 
presence of hydrocarbons. 

FIGS. 6 and 7 are flow charts of algorithms for estimation 
of concentrations of fluid components and for determination 
of presence of a particular fluid component, where a coated or 
doped SWCNT network is used. 

FIG. 8 graphically illustrates recovery time of a measured 
electrical parameter (here, current) to presence of ultraviolet 
light. 

FIG. 9 is a flow chart of a procedure for practicing the 
invention. 

FIGS. 10 and 11 schematically represent systems for prac- 
ticing the invention. 

FIGS. 12 and 13 graphically illustrate behavior of a mea- 
sured electrical parameter and response time with varying 
temperature T. 

FIG. 14 illustrates a model for estimation of location of 
upper and lower surfaces of a liquid, such as an oil-water 
mixture. 

FIG. 15 graphically illustrates behavior of a measured elec- 
trical parameter AR/R 0 for varying conc(N0 2 ) with varying 
relative humidity RH, fora SWCNT coated with hydroxypro - 
pyl cellulose. 

FIG. 16 graphically illustrates variation of an approxima- 
tion function for a measured electrical parameter with varying 
relative humidity RH for an SWCNT coated with an unspeci- 
fied substance. 

FIG. 17 schematically illustrates a downhole wireless 
transmission system for use with the invention. 

FIG. 18 is a flow chart of a procedure for practicing the 
invention for a selected environmental parameter. 

DESCRIPTION OF BEST MODES OF THE 
INVENTION 

FIG. 1 illustrates a nanosensor 11, constructed according to 
an embodiment of the invention, including 32 interdigitated 
electrode fingers (“IDEFs”), 12A-1 through 12D-8, con- 
nected to first and second voltage sources, with a controllable 
voltage difference. In FIG. 1, a 1 cmxl cm chip with Q 
sensing elements and SWCNTs bridging the gaps between 
adjacent IDEFs (Q=2-32). Each IDE is fabricated using con- 
ventional photolithography with a nominal finger width of 1 0 
\im and gap sizes d=8, 12, 25 and 50 pm, or any other desired 
set of gap sizes. The fingers and contact pads are thermally 
evaporated Ti (20 nm thickness) and Pt (200 nm thickness) on 
a layer of Si0 2 , thermally grown on a silicon substrate. The 
sensing material is bulk-produced SWCNTs from a HiPCo 
process (Rice University.), purified to remove amorphous 
carbon and metal impurities according to a procedure 
described in the literature. The purified SWCNTs are dis- 
persed in a selected solvent, such as dimethyl formamide, to 
create a suspension of SWCNTs. The suspension is soni- 
cated, then drop -deposited onto the interdigitated area of the 
electrodes. After the solvent evaporates, the SWCNTs forms 
a network connecting two adjacent IDEFs with a correspond- 
ing voltage difference. Any residue of the solvent was thor- 
oughly removed by drying the sensor under vacuum. The 
SWCNT density in the network can be varied by varying the 
amount of suspension or the concentration of SWCNTs in the 
suspension placed on the sensor area. This process produces 
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a statistically meaningful number of SWCNTs bridging the 
gap between two terminals to give reproducible performance. 
Alternatively, the SWCNTs can be grown directly on the 
electrodes using chemical vapor deposition. A selected elec- 
5 trical parameter response value, such as conductance or resis- 
tivity or electrical current or voltage difference, is provided or 
measured for the uncoated and coated SWCNT networks. 
Differences between the sensing elements are compensated 
for by varying a base resistance located in series with each 
to element. 

SWCNTs yield different signal responses when exposed to 
different fluids and vapors, and one must use pattern recog- 
nition or intelligent signal processing techniques for the iden- 
tification of the fluid constituent of interest. SWCNTs do not 
1 5 respond to exposure to certain fluids and vapors, and in those 
cases, coating or doping of the nanotubes may elicit a signal. 

In one version of the system shown in FIG. 1, eight of the 
IDEFs (12A-1 through 12A-8) are doped with mono-layer 
protected clusters (“MPCs”) ofAu, eight of the IDEFs (12B-1 
20 through 12B-8) are coated with hydroxypropyl cellulose, and 
sixteen of the IDEFs (12C-1 through 12D-8) are “bare” 
CNTs. In another version, each IDEF, or subsets including 
more than one IDEF, is coated or doped with a (different) 
selected coating or dopant material so that the presence (or 
25 absence) of up to 31 or 32 different target materials can be 
tested. A group of IDEFs, including a single IDEF, can be 
used to test for a given environmental parameter value (e.g., 
relative humidity, temperature, fluid pressure, fluid move- 
ment). Alternatively, a value for each of a selected group of 
30 environmental parameters can be determined, using an 
assembly of IDEFs illustrated in FIG. 1. 

In some initial tests, the electrical current (response value) 
through the sensors, at a constant voltage of 1 Volt, was 
monitored as different concentrations of chemicals, such as 
35 chlorine (Cl 2 ) and of hydrochloric acid (HC1) vapor were 
introduced to the sensor’s environment. A voltage difference 
of less than 1 Volt (or greater, if desired) can be used here. A 
computerized fluid blending and dilution system, Environics 
2040 (Environics, Inc. Tolland, Conn.), was used to create 
40 different concentration streams with a steady flow of 400 
cc/min during both exposure and purge periods. 

The electrical signal (current) was collected using a semi- 
conductor parameter analyzer HP41 55B (Agilent, Palo Alto, 
Calif.). Other equivalent electrical parameters, such as con- 
45 ductance or resistance, can be used as response values. In 
trials involving heating, a thermal controller, Micro -Infinity 
ICN77000 Series Controller (Newport Electronics, Inc., 
Santa Ana, Calif.) with a thermocouple, maintained a con- 
stant temperature for the sensor operation. Additionally, a 
50 vacuum pump and an ultraviolet lamp of wavelength 254 nm 
were employed on occasion to accelerate the recovery of the 
sensors between tests; other ultraviolet wavelengths, such as 
300 nm and 360 nm, can also be used to accelerate recovery, 
but are less effective than the 254 nm wavelength. 

55 Carbon nanotubes do-not sense presence of some fluids 
and vapors, due to the chemical and physical properties of 
CNTs as well as the nature of interaction between the fluid 
molecules and nanotubes. SWCNTs have been found to 
detect presence of NH 3 and/or N0 2 , based on the charge 
60 transfer between these fluids and SWCNTs. Early tests indi- 
cated that pristine SWCNTs do not respond at all when 
exposed to some industrial chemicals, such as chlorine, 
hydrogen chloride and methane. It is important to get some 
observable response before one can do signal processing or 
65 pattern recognition for selective identification. 

Carbon nanotubes coated with different polymers, such as 
chlorosulfonated polyethylene, hydroxypropyl cellulose, 
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polystyrene, polyvinylalcohol, etc. used in commercial poly- 
mer based chemical sensors available for organic vapor detec- 
tion, can provide specific interactions with a chemical species 
of interest. As this chemical treatment aims to provide a 
specific interaction between the carbon nanotube matrix and 
specific gas molecules, the treatment can improve the selec- 
tivity while maintaining the high sensitivity expected of a 
nano sensor. 

FIG. 2A is a schematic view of an IDE, used in the inven- 
tion, in which An assembly of SWCNTs, for example, is 
loaded with Pd (or another transition element, such as Pt, Ru, 
Rh, Ir, Os, Au, Ag and/or Hg) by sputter coating onto a pile of 
CNTs. The CNTs are then drop deposited onto the fingers in 
FIG. 2A and dried to create a sensor with CNT diameters of 
about 10 nm and an initial resistance in a range of 0.2-1 
kilo-Ohm. FIG. 2B is a scanning electron microscope (SEM) 
image in which the relatively bright spots are Pd particles 
(average diameter 1 0 mn) and the larger and darker bulk is an 
assembly of SWCNTs. A constant voltage of about 1 Volt is 
imposed across the CNT sensors, with CH 4 concentrations of 
6, 1 5, 30 or 100 ppm present to estimate the sensitivity of the 
IDEs in FIG. 2A to a hydrocarbon (here, CH 4 ) at room tem- 
perature. Test results exhibit an increase in current for 
increase in CH4 concentrations at these small values, with 
corresponding decrease in electrical resistance, or corre- 
sponding increase in conductance. A logarithm fit is found to 
most closely match the electrical properties data for most of 
the sensors for these small concentrations. 

FIG. 3 graphically illustrates sensor response (current in 
milli-Amps) versus time for each of several CH 4 concentra- 
tions (6, 20, 60 and 100 ppm), and also graphically illustrates 
fitting of a logarithmic function, 

V/V^a ln(0-6=ln{C7exp(6)}, (1) 

to relative response V/V 0 versus CH 4 concentration c, for 
different sensor networks. Resistance, as the electrical 
parameter, normally decreases with increasing concentration 
C of the constituent, while conductance, electrical current and 
voltage difference normally increase with increasing concen- 
tration c. The algorithm set forth in Appendix 1 and illustrated 
in Eqs. (A3)-(A6) can be applied to estimate concentration of 
one or more of the fluids CH 4 , C m H„ and/or CO x , by replacing 
the concentration C 1 or C 2 or C 3 =C by the quantity 

x=ln{C7exp (£>)}, (2) 

where the parameters a and b will vary with the particular gas 
constituent of interest. 

Methane, in the presence of the SWCNT network, may 
form a complex such as FI [Pd] .CH 3 . The H atoms in CH 4 tend 
to attract electrons from Pd, which in turn can obtain electrons 
from the SWCNTs to facilitate formation of the complex. 
This behavior should also be manifest for some or all of 
similar transition metals, such as Pt, Ru, Rh, Ir, Os, Au, Ag 
and/or Hg. The detection lower limit for CH 4 at room tem- 
perature, using a Pd-doped SWCNT network, is estimated to 
be a few hundred ppb to a few ppm. This compares with a CH 4 
detection lower limit Of 0.5-1 percent for conventional sen- 
sors, at temperatures T(min)^450° C. 

A coated or doped (“loaded”) SWCNT array can be used to 
estimate one or more values of a group of environmental 
parameter EP, including local relative humidity (RH), local 
temperature (T), local pressure (p), direction and velocity of 
fluid movement (represented by a current vector J) and/or 
presence/absence of a target molecule associated with an 
underground reservoir of a fluid mineral, such as oil, natural 
gas and the like. 
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A comparison experiment was conducted on sensors using 
pure, unloaded SWCNTs exposed to different gas analytes, 
with results shown in FIG. 3. The pure or unloaded SWCNT 
sensors showed no observable response to exposure when 
5 exposed to Cl 2 or to HC1, but displayed positive response 
signals, varying monotonically with concentration, for NO p , 
nitrotoluene, and benzene, and showed a negative monotonic 
response signal (opposite polarity) for NH 3 . Thus, the 
SWCNT sensor has some low level of discriminating power 
to by itself for exposure to some, but not all, gases. It is clear 
from FIG. 4 that SWCNTs sensors are sensitive to NO^, 
nitrotoluene, benzene and NH 3 , as demonstrated previously. 
The opposite polarities of sensor responses for NO^ and for 
NH 3 are believed to be due to the electron-acceptance and 
15 electron-donation characteristics of these fluids. In relative 
terms, SWCNT sensors have lower sensitivity to benzene and 
to nitrotoluene vapors than to NO^. 

A similar comparison experiment was carried out on the 
polymer-coated SWCNT sensors exposed to different ana- 
20 lytes, with results shown in FIG. 4. It is clear that addition of 
chlorosulfonated polyethylene coating to SWCNTs increases 
the sensitivity dramatically for exposure to 5 ppm Cl 2 com- 
pared to the pure SWCNT sensor, which shows no observable 
response even at 50 ppm of Cl 2 concentration. Although there 
25 is a signal for N0 2 from a chlorosulfonated polyethylene- 
coated SWCNT sensor, the sensing signal ratio of Cl 2 to N0 2 
is much larger than 1 . Therefore, this polymer-coated sensor 
can be used in a sensor array to give a distinct response signal 
for Cl 2 (and, similarly, for HC1) for discrimination between 
30 chemicals. 

These studies also show that both pure (unloaded) 
SWCNTs and chlorosulfonated polyethylene-coated 
SWCNTs do not respond to concentrations of up to 100 ppm 
HC1 fluid in nitrogen. Higher concentrations of HC1 were not 
35 tested as these high levels are not of interest for a nanosensor. 
In contrast, hydroxypropyl cellulose-coated SWCNTs 
respond to presence of HC1, but this sensor is also sensitive to 
N0 2 . Presence of the OH groups in the polymer coating may 
be responsible for the response signal differences in interac- 
40 tion with acidic fluids. Because this sensor gives a significant 
response to HC1 where other SWCNT sensors do not, the 
hydroxypropyl cellulose sensor can be used in a sensor array 
to provide a chemical signature that differentiates the HC1 
fluid from other chemicals. 

45 We have demonstrated a simple nano -chemical sensor 
using polymer-coated SWCNTs as the sensing medium. 
Because pristine (unloaded) nanotubes do not respond 
observably to some chemicals of interest, it is important to 
explore coating or doping techniques to promote observable 
50 responses so that a broad application coverage can be 
ensured. We have found that the polymer coating enables 
selective sensing of, and discrimination between, chlorine 
and hydrochloric acid vapor at a sensitivity level of 5 ppm and 
above. It is important to recognize that coating or doping 
55 alone is unlikely to provide absolute discrimination. It is 
recommended that the two types of loadings be used together. 
As with most sensors of any size or exploiting any property 
change, pattern recognition techniques are a valuable and 
necessary complement to provide discrimination. In that 
60 regard, the use of sensor arrays with multiple elements is an 
effective approach to chemical sensing, wherein the data from 
multiple sensors can be routed to a signal processing chip, 
integrated into the system, for data fusion and analysis. 
Advanced signal processing and pattern recognition tech- 
65 niques can be used to confirm (or refute) the assumed pres- 
ence of a given species, in addition to the help from the 
selective coatings. Sequential use of multiple sensing element 
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arrays offers additional operational freedom, when sensor 
recovery is slow and is a rate limiting process. Under such 
circumstances, a sensor would always be available while 
other sensors are in recovery mode. 

Using the results shown in FIG. 4, a procedure, shown as a 
flow chart in FIG. 5, can be implemented to detect presence of 
one, two or more target molecules, such as halogens, C m H w 
and/or an oxygen-containing compound (alcohol, ketone, 
aldehyde, etc.) C n f{ p O q in an unknown fluid G (which may 
contain N0 2 , .NH 3 and one or more of the target fluid mol- 
ecules). In step 51, a known increment of a selected (possible) 
component of the fluid G (reference component, such as N0 2 
or NH 3 , or target component, such as CH 4 or C 2 H 6 ) is added 
to and mixed with a portion of the fluid G to provide an 
augmented fluid G\ In step 52, response values, V 0 (G;meas) 
and V 0 (G';meas); are measured for the uncoated SWCNT 
sensor In step 53, the SWCNT network is loaded (coated or 
doped) with a selected coating or dopant no. nl (e.g., Pd, Pt, 
Rh, Ir, Ru, Os, Au, Ag and/or Hg), is exposed to the fluid G 
and (separately) to the fluid G', and the corresponding 
response values V Ml (G;meas) and V wl (G';meas) are measured 
or otherwise provided. In step 54, the SWCNT sensor is 
loaded (instead) with a different selected coating or dopant 
no. n2, is exposed to the fluid G and (separately) to the fluid 
G', and the corresponding response values, V„ 2 (G;meas) and 
V W 2 (G';meas) are measured or otherwise provided. In step 55, 
the system uses a selected algorithm to estimate the concen- 
tration (e.g., in ppb), or concentration upper bound or con- 
centration lower bound, of one or more of the components 
believed to be present in the fluid G (or in the fluid G’). Two 
suitable algorithms for determination of fluid component 
concentration are disclosed in flow charts in FIGS. 6 and 7, 
discussed in Appendices A and B, respectively. 

It is assumed initially in Appendix A that (i) the response 
value difference varies substantially linearly, over a suitable 
concentration range, with concentration difference of a single 
constituent that is present and (ii) the response value differ- 
ence, in the presence of two or more fluid constituents in the 
fluid G, is the sum of the response value differences of the 
single constituent fluids. Linear response coefficients a.^ for 
the response value differences are assumed to be determined 
experimentally or otherwise provided. As an example, 
assume that one reference fluid (e.g., N0 2 or NH 3 ) plus first 
and second target fluids (e.g., CH 4 and C 2 H 6 ), are suspected 
to be present in the fluid G. Estimates of each of the concen- 
tration values c m0 for the initial (unaugmented) fluid G are 
obtained from inversion of an MxM matrix equation relating 
these concentration values to response value differences for N 
coatings, where M (=2) is the number of fluid components 
(reference and target) believed to be present and N (=1) is the 
number of loadings (coatings or dopants) used for the mea- 
surements. 

The approach discussed in Appendix A, and illustrated in 
FIG. 6, allows separate weights, w w and w' w (or the same 
weight), to be assigned to the measurements of the initial fluid 
and augmented fluid. Preferably, at least two of the weight 
values in Eq. (3) are positive (e.g., (w l5 w 2 ) or (w' l5 w' 2 ) or 
(w ls w' 2 ) or (w' l3 w 2 )) for the example with N=2, and the 
relative sizes of the non-zero weights reflect the relative 
importance of the response measurements. If, as is likely, the 
four response measurements are believed to be equally 
important, one can choose w 1 =w 2 =w' 1 =w , 2 =l. One can 
ignore one or two of the four measurements, in which event 
the corresponding weight value(s) is set equal to 0. 

The response coefficients abused in Eqs. (A1 ) and (A2) are 
not necessarily positive. For example, the response coeffi- 
cient a for the fluid constituent CH 4 is positive for several of 
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the SWCNT coatings used, while the response coefficient for 
NH 3 is observed to be negative for some of these coatings. 

Appendix B, and the corresponding flow chart in FIG. 7, 
discuss an approach for estimating an upper bound, or a lower 
5 bound, of a concentration value of a selected molecule (e.g., 
C m H w and/or C^H^O^), again assuming that a response value 
difference varies substantially linearly with a concentration 
value difference of the selected molecule. This bound is com- 
puted separately for each selected molecule. Appendices A 
to and B apply to loaded (coated or doped CNT networks. 

Exposure of the coated SWCNT network to ultraviolet 
light can reduce the recovery time (normally ten hours or 
more) required to return the network to a substantially 
uncoated condition, by promoting accelerated detachment of 
15 the coating material from the SWCNT network. FIG. 8 
graphically illustrates measured response to absorption of 
CH 4 versus time with exposure to ultraviolet light for two 
different amounts of the same coating (sensors 1 and 2), 
indicating the improvement in recovery time where ultravio- 
20 let light is applied to lower the energy barrier to desorption of 
the adsorbed chemical from the SWCNT network. The 
chemical and environmental time histories may cause the 
response curves for the same sensor to differ. 

For some relatively small molecules, such as methane 
25 (CH 4 ), ethane (C 4 H 6 ), other hydrocarbons, and oxides of 
carbon (CO x ; x=l-2), an SWCNT network, doped with a 
transition element (“TE”,) such as Pd, Pt, Rh, Ir, Ru, Os, Au, 
Ag and Hg, is used to detect presence of these molecules by 
detecting a change in an electrical parameter (conductance, 
30 resistance, current or voltage difference) or response value 
associated with a path defined by an SWCNT network, which 
extends between two electrodes having a controllable voltage 
difference or current. Some molecules, including nitrotolu- 
ene and phenol, are relatively strong electron donors and/or 
35 electron acceptors, and these molecules’ presence can be 
readily detected using “bare” or unmodified SWCNTs. Other 
molecules, including but not limited to methane, other hydro- 
carbons and carbon oxides, manifest little or no electron 
donor or electron acceptor action so that monitoring an elec- 
40 trical parameter value V of an unmodified or “bare” SWCNT 
network will not, by itself, indicate presence or absence of 
these molecules. 

Where CH 4 is adsorbed in a SWCNT/Pd matrix, the com- 
bination forms a weakly bound complex, such as Pd 6 (CH 4 ) -6 , 
45 where 6 is a relatively small positive number that need not be 
an integer. Methane, other hydrocarbons and carbon oxides 
are “greenhouse” gases and require detection capabilities in 
the ppb-ppm range to have much utility in environmental 
monitoring. The sensing platform is similar to that illustrated 
50 in FIG. 1, where the SWCNT/TE compound serves as a 
current or voltage modifier whose electrical parameter value 
V changes when molecules of a selected target chemical, such 
as CH 4 or C m H„ or CO x , are adsorbed on the SWCNT/TE 
compound, and the magnitude of the value V is monotoni- 
55 cally increasing (not necessarily linearly) with increase in the 
amount of the target chemical present. This approach for 
detecting presence of a target molecule is often at least one 
order of magnitude more sensitive at room temperature than 
is detection using catalytic beads or metal oxides, in part due 
60 to (1) nano scale-induced charge transfer between TE-loaded 
nanotubes and CH 4 molecules and (2) a relatively large sur- 
face area per unit volume for SWCNTs, which permits 
enhanced adsorption of hydrocarbon molecules (e.g., CH 4 
and C m HJ. 

65 Fabrication of a sensing platform for the SWCNT/TE net- 
work begins with sputter coating of about 10 nm thick Pd (or 
another suitable TE) onto a pile of SWCNT powder. The 
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TE-loaded SWCNTs are then dispersed into distilled, deion- 
ized water (e.g., 0.1 mg of SWCNT/Pd in 1 0 ml of the water). 
This solution is then sonicated and drop deposited onto inter- 
digitated electrode fingers to create an electrical sensor with 
an initial resistance in a range of about 0.2-1 kilo-Ohm. Cur- 
rent through the network, with a 1 Volt difference, was moni- 
tored where 6,15,30 and 1 00 ppm of CH 4 was present, using 
a fluid stream flow of about 400 cc/min, during exposure and 
during purge. Apart from preparation of the doped or loaded 
SWCNT network (as distinguished from coating a SWCNT 
network), the procedure for estimation of constituents present 
in a fluid using a doped SWCNT network, is parallel to the 
procedure using a coated SWCNT network set forth in the 
FIG. 5 flow chart. 

Vacuum pumping and exposure of the network to ultravio- 
let light (X=254 nm) are used to reduce the recovery time of 
the SWCNT/Pd (i.e., removal of the CH 4 ) between tests (no 
recovery, if these recovery acceleration procedures are not 
implemented). FIG. 8 graphically illustrates measured 
response for CH 4 versus time with ultraviolet light present, 
indicating the improvement in recovery time where ultravio- 
let light is applied to lower the energy barrier to desorption of 
CH 4 molecules from the SWCNT/Pd network. 

FIG. 9 is a flow chart illustrating a procedure for practicing 
an embodiment of this invention. In step 91, an array of one or 
more SWCNTs is loaded (coated or doped) with one or more 
selected chemical substances, chosen so that an electrical 
parameter value (ELPV), measured along at least one loaded 
SWCNT, varies monotonically with variation of a value of a 
corresponding selected environmental parameter EP. In step 
92 (optional), the value of the EP is allowed to vary. In step 93, 
a measured electrical parameter value ELPV is provided cor- 
responding to the value of the EP. In step 94, the EP value, 
corresponding to the ELPV provided, is measured or esti- 
mated. In step 95 (optional), the ELPV and/or corresponding 
EP value is transmitted. 

The EP may be one or more of relative humidity, tempera- 
ture, fluid pressure, fluid level (for a liquid), direction and rate 
of fluid flow (represented by a fluid current vector J) and/or 
presence/absence of one or more target molecules in the fluid. 
FIG. 10 schematically illustrates a chip system 101 that can 
be used to practice the invention. The apparatus includes one 
or more of: a chip 102 configured with one or more nano- 
sensors to sense the presence/absence of one or more target 
molecules; a relative humidity sensor 103 to sense local 
humidity; a temperature sensor 104 to sense local tempera- 
ture; a pressure sensor 105 to sense local pressure; an acoustic 
sensor 106 to sense fluid level (for a liquid) adjacent to the 
sensor; a fluid flow sensor 107 to sense direction and rate of 
fluid flow (represented by a current vector J); and a wireless 
(or hard wired) transmitter or transceiver 108 (optional) to 
transmit sensor data to a database. One or more of the systems 
101 is positioned in each producing horizon of a fluid mineral. 

Where the mineral is oil or another mineral -containing 
liquid, steam or water or another inert liquid substance can be 
injected at a mineral -producing horizon (a vertical or depth 
range within which the mineral is located, isolated from any 
adjacent producing horizon) to move or displace the mineral 
and/or to provide additional pressure to bring the mineral to 
the surface for recovery. A fluid flow sensor can sense the 
direction and velocity of movement of the mineral, for pur- 
poses of monitoring the effectiveness of the steam pressur- 
ization. If, for example, the liquid mineral is not moving in an 
appropriate direction, the origin of the steam pressurization 
can be changed. 

FIG. 11 schematically illustrates injection of steam or 
another relatively inert fluid into a reservoir R, or a suspected 
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reservoir, to estimate relevant physical and chemical content 
of the reservoir and movement of a mineral -containing fluid 
in response to this injection. One or more sensor chips 102 
(FIG. 10) is positioned in a mineral producing horizon H to 
5 sense one or more of chemical composition, relative humid- 
ity, fluid temperature, fluid pressure, fluid flow direction, and 
estimated amount or concentration of mineral present adja- 
cent to the chip 102. A transceiver 108, positioned on or 
associated with the chip 102 in FIG. 10, transmits signals 
to representing measurements made by the sensors. Preferably, 
a network of such chips provides sensor information that can 
be correlated between adjacent chips so that, for example, a 
fluid flow or velocity pattern can be estimated for relevant 
parts of the reservoir R. This flow pattern can be estimated for 
15 horizontal movement (only) of the fluid or for combined 
horizontal and vertical movement. 

FIGS. 12 and 13 graphically illustrate the environmental 
effects of varying temperature upon relative response V/V 0 
and upon response time, respectively, for SWCNTs plus Pd. 
20 The relative response for CH 4 in FIG. 12 appears to reach a 
relative minimum at a temperature T=T(min)«40° C. and to 
increase monotonically as IT-T(min)l increases. At room 
temperature, T®20° C., the relative response is approximately 
three times the relative response at T=T(min). The response 
25 time (required for approximate equilibration of the response 
after exposure to a fluid containing the target molecule) 
decreases monotonically with increasing temperature T and 
increases with the amount of load (coating or dopant) used to 
treat the SWCNT. 

30 Temperature can be measured using a temperature 
micosensor, such as a SIP Type, TD Series, available from 
Honeywell. This sensor should accurately measure, and pref- 
erably record, temperatures in a range from -40° C. up to 
150° C. 

35 Gas pressure can be measured using a pressure micro sen- 
sor such as a 140 Series Silicon Pressure Sensor, available 
from Honeywell. This sensor should accurately measure fluid 
pressures in a range from 0 psi to 1 50 psi. 

The S WCNT/TE sensor, with TE=Pd, has been tested at 1 5 
40 ppm and 30 ppm concentrations of CH 4 at temperatures of 
T=40° C., 80° C. and 150° C. The response parameter value 
(e.g., conductance or current) increases with increasing tem- 
perature, perhaps due to an enhanced catalytic effect of 
SWCNT/Pd binding with increasing temperature, as illus- 
45 trated in FIG. 13 for CH 4 . 

Experimental results for other hydrocarbons, for CO x , for 
alcohols, for ketones and for aldehydes are qualitatively simi- 
lar to those for CH 4 . 

U.S. Pat. Nos. 6,672,163, 6,817,229, and 6,938,458, issued 
50 to Wei Han et al, disclose a system for in situ characterization 
of a down-hole fluid, using a comparison of transmitted and 
received signals that are propagated through the fluid, which 
may contain drilling mud and other solids, water, oil and other 
liquids that need to be characterized. The method uses vary- 
55 ing frequency signals to estimate attenuation, scattering and 
other effects produced in a non-homogeneous liquid-solid 
environment. Below 22 percent (by volume) solid concentra- 
tion, attenuation a varies with frequency approximately as 
a=Af+B, where A and B are substantially independent of 
60 frequency f. In a pure liquid, the attenuation a=A , f t +B', with 
k«2. As the solid fraction increases, the exponent k decreases 
toward k=l. In a pure liquid, the signal propagates with 
attenuation a per unit length, is reflected by a solid, unmoving 
barrier, and returns to a receiver. 

65 Consider a model of liquid (e.g., oil) level estimation illus- 
trated in FIG. 14, where a source-sensor 141 is located at a 
height hi above an exposed (upper) surface SI of a liquid 142 
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having a depth h2. An acoustic or electromagnetic signal 142 
is provided by the source-sensor 141 and propagates in an 
ambient medium 143 to the upper surface SI. A calculable 
first reflection coefficient pi (J) of the signal energy is 
reflected at the surface SI and returns to, and is measured at, 5 
the source-sensor 141, with a total attenuation, pl(/)*exp{- 
2 al(/)-hl}, that is characteristic of signal propagation by a 
distance 2/zl. Here, ol(f) is a frequency-dependent attenua- 
tion coefficient for the ambient medium 142. By measuring 
the attenuation at two distinct frequencies, f=fl and f=f 2 , and 10 
using known attenuation coefficients, al(fl) and al(f 2 ), one 
can estimate the present height hi of the source-sensor 141 
above the upper surface SI of the liquid 144. By appropriately 
gating in time of the initial signal produced and the return 
signal from the upper surface 51 , the incremental signal from 1 5 
a portion of the propagating signal that has passed through the 
liquid 144 can be deleted. 

A second portion (calculable refraction coefficient p2 (/)) 
of the initial signal passes through the upper surface SI, 
propagates an unknown distance li2 through the liquid 144 20 
with a known attenuation coefficient a 2 (/), and is reflected at 
a lower surface S2 (assumed to be planar) of the liquid 144 
with a calculable reflection coefficient p3 (/), and propagates a 
distance h2 to the upper surface SI of the liquid 144. A fourth 
portion (calculable refraction coefficient p4 (/)) of the signal 25 
passes through the upper surface SI and propagates a distance 
hi to the source-sensor 141. From well known electromag- 
netic signal relations at an interface between two materials at 
approximately 0 ° incidence angle, the quantities pi (/), p 2 (/), 
p3 (/) and p4 if) are estimated to be 30 


P 1 (J)={(n2(f)-n(J))l(n2(J) + n\ (f))} 2 , 

(3A) 

P 2(f)=l- P m 

(3B) 

p3(f)={(n3(f)-n2(f))/(n3(f)+n2(f))} 2 , 

(3C) 

pw=pm. 

(3D) 


where nl (/), n 2 (/) and n3 (J) are the refractive indices of the 
ambient medium (wellhead space), liquid 152 and lower sur- 40 
face (S2) material, respectively, which are known or deter- 
minable. 

The attenuated portions of the two signal components thus 
become 

sl=pl-exp{-2al(/)-/zl}, (4) ^ 

for the portion reflected (once) at the upper surface SI, and 

s2=p2p3p4-exp{-2al(f)-hl-2a2(f)-h2}+sl (5) 

50 

for the portion that passes through the upper surface SI and 
through the liquid (once in each direction). The distances or 
depths, hi and h2, can be determined from Eqs. (4) and (5). 

The information disclosed in the Wei Han et al patents can 
be adapted to provide an estimate of distance hi to an upper 55 
surface of the liquid and a representative depth h 2 of the liquid 
for the producing horizon. 

Location of the upper surface SI in FIG. 14 can also be 
estimated as follows. Two or more adjacent nanosensors, 
attached to or associated with a drill bit or other instrument 60 
whose elevation is controllable and measurable, are lowered 
into the head space of the well site. An electrical property, 
such as conductance, is measured at a sequence of times. 
When the two nano sensors are in the headspace, the conduc- 
tance has, a first (known) value. When the two nano sensors 65 
move past the first surface SI and into the liquid 144 in FIG. 

14, the conductance abruptly changes to a second value that is 
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distinct from the first value. The upper surface SI is located at 
the vertical height at which this abrupt change occurs. Loca- 
tion of the lower surface S2 would require other measure- 
ments and computations. 

Effects of varying relative humidity (RH) on relative 
response, referenced to the response (V 0 ) at RH=0 percent, 
have been measured for several coating materials and doping 
materials, for a sequence of RH values. It is expected that 
relative response V/V 0 will decrease monotonically as the RH 
value increases, in part because the presence of a polar sub- 
stance such as water would interfere with, and partly mask, 
the change AW in an electrical parameter, for substantially all 
coating and dopant materials of interest. This expectation is 
borne out in measurements of AR/R 0 , at T=T 0 =23° C., pre- 
sented graphically in FIG. 15, for varying conc(N0 2 ) with the 
coating material hydroxypropyl cellulose, applied to a 
SWCNT, for RH=0, 15, 30, 50, 70 and 90 percent. 

More generally, the curves for relative response V/V 0 can 
be approximated or estimated by a parametrized curve (FIG. 
16), such as 

V/V 0 (RH; T 0 )=F e ((RH/RH 0 ) m ;a)~a ■ seek { (RH/ 

RH 0 r}+(l-a) ( 6 ) 

where RH 0 is a reference RH value, m is a positive number 
and 0 <a^l, each value being chosen for the particular load- 
ing material (coating or dopant) of interest. For sufficiently 
small values of the quantity RH/RH 0 , the parametrized curve 
in Eq. ( 6 ) is further approximated as 

V/V 0 (RH; T 0 )~'a { 1 - (RH/RH 0 ) 2m /2}+(l-a), 

=l-(a/2)(RH/RH 0 f m , (7) 

which is linear and decreasing in the quantity (RH/RH 0 ) 2m , 
and thus linear in the variable RH if 2 m=l . More generally, 
the measured relative response value V/V 0 (RH;T 0 ) is sub- 
stantially monotonically decreasing in the value RH and 
resembles a trapezoid with a non-zero tail value, as illustrated 
in the generalized curve shown in FIG. 16. 

Introducing the dimensionless variable 

u=(RH/RH 0 ) m , (8) 

several interesting values of u can be identified on FIG. 1 6: ( 1 ) 
a knee value, u=w ia , where the first w-derivative is changing 
most rapidly; ( 2 ) a maximum slope value, u=u^, where the 
(negative) slope or u-derivative is a maximum value; and ( 2 ) 
a second knee value, u=u ^ 2 (>u fl ), where the first u-deriva- 
tive is again changing most rapidly. The values are estimated 
from the identities 

(3 2 F c /3« 2 ) wA: =a -sech 3 u{ sin h 2 u- 1 }^=0, (9A) 

(d 3 F y3« 3 ) wMjS -=a-5 , ecfc 4 «{cos k 2 u-4} MS =0, (9B) 

The value of the parameter a in the interval 0<a^l is not 
relevant here. The solutions of Eqs. (9 A) and (9B) are verified 
to be 


u K =smh 1 {l)=ln{V2±l }, 

(10A) 

w xw5 = l n (2±'3}, 

(10B) 


With appropriate choices of the parameter values RH 0 and m 
in Eq. ( 8 ), the experimentally observed locations of the values 
u ms an d u ^ can be matched. 
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An appropriate value of the parameter a can be determined 
as follows. Let u=u^ correspond to the value of u for which 
F e (u / )=g (0<g<l). From Eq. (6), this requires that 


a -seek (uj )+ 1 -a =g, 

(UA) 

sech(uy )= 1 - ( 1 -g)/a), 

(UB) 


which requires that a lie in a reduced range, l-g<a^i 1 . Equa- 
tions (9A), (6B) and either (10A) or (10B) can be used to 
estimate RH 0 , m and a. 

If the approximation in Eq. (6) is adopted, the zero point 
relative response V/V o (RFl=0) for a particular coating or 
dopant material can be compensated for the presence of mois- 
ture (RH>0) by a compensation factor such as 

V o {RH=())=V{RH>0)l{a-sech{{RH/RH o ) m +{l-a)}. (12) 

Other approximations, such as replacing the sech(u) function 
in Eq. (6) by sec(u), or by exp {-u 2 /u 0 2 }, for example, can also 
be used here for RH compensation. 

For a given coating or doping material and fixed tempera- 
ture, the measured relative response V(RH;meas)/V 0 , as rela- 
tive humidity RH is increased over a sequence of values, can 
be compared with corresponding reference values V(RH;ref; 
h)/V 0 for each of a plurality of candidate fluids (h=l, . . . , H) 
to determine if a particular candidate fluid is present. Appen- 
dix C sets forth an analytical procedure for determining if a 
target fluid is likely to be present, from a comparison of 
measured relative response values V/V 0 for a variable envi- 
ronmental parameter, such as relative humidity, temperature 
or pressure. 

This technology combines physical and chemical measure- 
ments in a single system and allows simultaneous measure- 
ments to be made for real time monitoring. By providing a 
system, shown schematically in FIG. 10, in each producing 
horizon of a plurality of adjacent mineral wellsites, real time 
monitoring of an entire mineral field, and of its reaction to 
pressurization or other perturbation at one or more locations, 
can be implemented; real time feedback is now available for 
modification of the perturbing force(s) in the field. 

A nano sensor array has been developed to estimate a ratio 
of H 2 0 (vapor or liquid) to a representative hydrocarbon(s) 
(“HC ,” for example, a paraffin) in a wellsite fluid, in or adja- 
cent to a headspace above a drill bit or other exploratory 
mechanism. For a specified electrical parameter ELP, such as 
change in electrical resistance, in electrical conductance, in 
voltage difference or in electrical current, a sequence of ref- 
erence ELP values (AV=ELPV) is estimated or otherwise 
determined, as a function, AV^FljVP^H-jO)}, of H 2 0 
vapor pressure present and, separately, as a function, 
AV 2 =F2{PV(T;HC)}, of HC vapor pressure present in or 
adjacent to the headspace. The vapor pressure VP of each of 
the H 2 0 and the HC mixture will increase monotonically with 
increasing temperature T. Water molecules are highly polar, 
whereas an HC molecule (e.g., C 6 H 6 ) is generally nonpolar. 

For a given wellsite and associated well headspace, one or 
more measurements are performed of the ELPV AV, with 
contributions to AV from H 2 0, from drilling mud (“DM”), 
and from the extant HC mixture. For a given (measured) 
wellsite temperature T, the measured total vapor pressure 
VP(T) in the headspace is a sum of partial fractions, 

VP(D=aVP(7:-/f 2 0)+pVP(2:-/fC)+(l-a-P)VP(2:- 

DM) (13) 
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of the component fluids, where VP(T;DM) is often small 
enough to be ignored, equivalent to setting l-a-|3~0. The 
ratio 

A F/A V 2 =F1 {VP (T;H 2 0)}IF2 {PV(T;HC) }~R (T)-al 
5 (l-a) (14) 

is taken as a measure of the H 2 0/HC ratio present in the head 
space, where R(T) is known or determinable. The two quan- 
tities AV X and AV 2 are both positive, which ensures that a lies 
to in the range 0<a^l. The quantity VP(T;H 2 0) can be esti- 
mated from empirical or theoretical data, and total VP can be 
measured so that VP(T;H 2 0) and VP((T;HC)=VP(total)-VP 
(T;H 2 0) can be estimated for the arguments of the functions 
FI and F2. 

1 5 A wireless transmission system for data accumulated in a 
down-hole environment is disclosed by Jing Li in “A Wireless 
Portable Carbon Nanotube-Based Chemical Sensor System,” 
unpublished. FIG. 17 schematically illustrates components of 
the system 170 , which includes a chemical nanosensor mod- 
20 ule 171 , whose output signals (four, in one embodiment) are 
received by a multiplexer 172 . The MUX output signal is 
augmented by a constant current source 173 and is passed 
through a buffer amplifier 174 (optional) and received by a 
microcontroller 175 (e.g., a PIC16F688, available from 
25 Microchip Technologies). The microcontroller 175 also 
receives a signal from at least one of a temperature micro sen- 
sor 176 , a vapor pressure microsensor 177 , a flow microsen- 
sor (direction, velocity) 178 , a relative humidity micro sensor 
179 and an acoustic or electromagnetic distance sensor 180 . 
30 The microcontroller output signals (e.g., in RS-232 format) 
are received by a wireless server 181 (e.g., a MOXA W2150), 
which transmits these signals, using 802.11b protocol, 
through an antenna 182 to a LAN. The transmitted signals are 
received by a receiver (not shown) that is spaced apart from 
35 the transmission system 171 , for display and/or further signal 
processing. 

A nanosensor includes two or more functionalized CNTs 
for sensing presence, and optionally amount, of one or more 
target chemicals. A microsensor may use nanosensors and/or 
40 conventional sensors to estimate one or more physical param- 
eters, such as temperature, fluid pressure, relative humidity 
and distance. 

Appendix A. Estimation of Concentration of a Fluid Compo- 
45 nent. 

FIG. 6 is a flow chart of a procedure for estimating con- 
centration of a fluid component that is suspected to be present. 
This approach requires knowledge of certain response coef- 
ficients relating change of a response value to increase or 
50 decrease in concentration of a particular fluid component G. 
In step 61 , response coefficients a w m (n=l, . . . , N; m=l, 
2, . . . , M) are estimated or otherwise provided for M single 
constituent fluids (numbered m=l, . . . , M; M^2) (such as 
N0 2 , Cl 2 and HC1 or for any other three fluids of interest), in 
55 the presence of SWCNTs individually coated with each of N 
selected coatings (numbered n=l, . . . , N; N^l), for example, 
with chlorosulfonated polyethylene (n=l) and (separately) 
with hydroxypropyl cellulose (n=2). 

For an SWCNT coated with coating number n, 

F„(G;meas)- F 0 (G;meas)=a„ ?1 m c N02 +a n2 c cl2 +a n ? 

3 c HCb (Al) 

where, for example, c NQ2 represents the N0 2 concentration 
(e.g., expressed in ppm or in ppb). The set of response coef- 
65 ficients {a w ^} m , for different coatings, n=nl and for n=n2 
(*nl), will differ from each other, but each set is determined 
or estimated by measurement of the response value differ- 
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ence, V w (G;meas)-V 0 (G;meas), of coated (n) versus 
uncoated (n=0), for each of the two (or, more generally, 
Mi^2) single constituent fluids present in a known concen- 
tration. For a single constituent fluid N0 2 and no coating 
(n=0), for example, a o l =0.034±0.002. 5 

In step 62, a known increment of one (or more) of the 
(suspected) constituent (e.g., N0 2 or Cl 2 or HC1), is added to 
the unknown fluid G to provide an augmented fluid G\ In step 
63, the response values, V 0 (G;meas) and V 0 (G';meas), for the 
uncoated S WCNT network (n=0), in the presence of the fluids 1 0 
G and G', are measured or otherwise provided. In step 64, the 
response values, V w (G;meas) and V w (G';meas), for the 
S WCNT network coated with the (single) coating number n, 
in the presence of the fluids G and G', respectively, are mea- 
sured or otherwise provided. 15 

In step 65, an error function €, defined by 

2£(x,y,z)=1. n w n -{ F„(G;meas)- V 0 (G;meas)-a nl c l -a ri , 

C2+«„ 5 3C 3 } 2 +2„w ■{ U„(G’;meas)- Uo(G';meas)- 

iCi+Aq)-fl w ^c 2 +fl nj 3C 3 } 2 , (A2) 20 

is provided, where c l5 c 2 and c 3 refer to the concentrations of 
the reference molecule, the first fluid molecule and the second 
fluid molecule, Ac x is a known concentration increment of a 
selected one (cj of the reference molecule, the first fluid 2 5 
molecule or the second fluid molecule, added to the fluid G to 
provide the augmented fluid G', and w w , and w' w are selected 
non-negative weight values. The two sums in Eq. (A2) rep- 
resent the contributions of the initial composition and the 
augmented composition, respectively. These sums over nmay 30 
include one, two, three or more loadings for which the 
response coefficients are known. In this example, n=l, 2. 

The error function €(c x , c 2 , c 3 ) is to be minimized with 
respect to choices of the (unknown) concentration values c l5 
c 2 and c 3 . Differentiating € with respect to each of the vari- 35 
ables c : , c 2 and c 3 , in step 66, one obtains three coupled linear 
equations in these variables: 


{(w : + wi )aj A + (w 2 + W 2 )al l }c l + 

{(w! + wi)ai 4 fli j2 + (w 2 + W 2 )a 2 ,ia 2 ,2}c2 + 

{(wi + wi)a U fli >3 + (w 2 + W 2 )a 2 ^a 2 } 3 }c 3 = 

M (V[ ( G\meas ) - Vq (G;meas)) + 

>^1(^1 (G'\meas) - Vq{G' \ meas)) - 0 i,iAcx)}ax,i + 
{w 2 V 2 (G;meas) - V§{G\meas)) + 

W 2 {V 2 {G' \me as) - Vo(G';meas)) - 
02,1 Aci )}a 2rl , {(wx + W i )a lA a U2 + (w 2 + W 2 ) 

02,102,2^1 + {(Wx + w[)a 2 l2 + (w 2 + w 2 )a\ 2 }c 2 


(A3) 40 


45 


50 


In step 66, Eqs. (A1 )-(A5) in the unknowns c 1 , c 2 and c 3 are 55 
determined, using standard matrix inversion techniques, after 
verification that a 3x3 (more generally, MxM) coefficient 
matrix for the vector [c l c 2 c 3 ]* r has a non-zero determinant. 
These solutions, [c 1 c 2 c 3 ] tr , provide estimates of the concen- 
tration values of the corresponding chemicals in the fluid G 60 
(or in the fluid G’) in step 66 of FIG. 6. 

Preferably, at least two of the weight values in Eq. (A2) are 
positive (e.g., (w l5 w 2 ) or (w' 1? w' 2 ) or (w 1; w' 2 ) or (w’ l5 w 2 )), 
and the relative sizes of the non-zero weights reflect the 
relative importance of the response measurements. If, as is 65 
likely, the four response measurements are believed to be 
equally important, one can choose w 1 =w 2 =w , 1 =w , 2 =l. One 


can ignore one or two of the four measurements, in which 
event the corresponding weight value(s) is set equal to 0. 

The response coefficients a w m used in Eqs. (Al) and (A2) 
are not necessarily positive. For example, the response coef- 
ficient a* . for the fluid constituent N0 2 is positive for several 
of the SWCNT coatings used, while the response coefficient 
a y for NH 3 is observed to be negative for at least one of these 
coatings. 

More generally, where M reference fluid components 
(numbered m=l, . . . , Ml) and target fluid components 
(numbered m=Ml+l, . . . , M1+M2=M) with unknown con- 
centrations are believed to be present and N coatings (num- 
bered n=l, . . . , N), the error function s (analogous to Eq. 
(A2)) is defined by 


2e(cx, . . . , c M ))=p„w fJ -{F P fJ (G;meas)-F 0 (G;meas)- 
?i m a n ^c m } 2 +l. n w' n -{ F„(G';meas)- Fo(G’;meas)- 
2 ra 0„,x(c ra +Ac m )} 2 , (A6) 

where one, or more then one, concentration value c m is aug- 
mented by a known amount Ac m . The error function € is 
minimized by differentiation with respect to each of the 
unknown concentration values c m . This yields M coupled 
equations 


2 w w w T v «( G ; meas )- v o( G ; meas )- 2 m a «^ c ^o} a «^o+ 

2 w w , „-{V„(G , ;meas)-Vo(G , ;meas)-2 ra a„^(c ra0 + 

Ac„, o )K,™o=0> (A7) 

for index values m0=l, 2, . . . , M. These can be restated in a 
matrix format as 

VoiGimeas^+^w^i^a^Ac^a^o, 

+Z„w'„-{ r„(G’;meas)- Fo(G’;meas)}. (A8) 

After verifying that the determinant of the MxM matrix of 
coefficients for the quantities c m0 in Eq. (A8) is non-zero, this 
MxM matrix can be inverted to determine estimates for the 
concentration values c m0 (m0=l, . . . , M). These concentra- 
tion value estimates will depend, in part, upon the relative 
values chosen for the weight values w„ and w' w for the coat- 
ings. Where one or more of the reference molecule concen- 
tration values c m0 (m0=l, . . . , Ml) are known in advance, the 
estimates for these reference concentration values can be 
compared with the corresponding known values to evaluate 
the likely accuracy of the remaining estimated values. 

The approach set forth in this Appendix A can also be used 
to estimate an initial concentration value c m0 where the CNT 
network is doped or otherwise loaded, rather than being 
coated. 

Appendix B. Determination of Bound on Fluid Component 
Concentration. 

A second algorithm does not require provision of a large 
number of response coefficients a zV but only seeks to deter- 
mine if a particular target molecule is present in at least a 
selected concentration. For a selected coating, such as chlo- 
rosulfonated polyethylene or hydroxypropyl cellulose or a 
transition element TE, on the CNT, a measurement of the 
response value difference AV=V(coated)-V(uncoated) is 
taken for modified fluids, G'(l) and G’(2), where each of two 
distinct supplemental concentration values, A^mO) and A 2 c 
(mO), respectively, for a selected molecule no. mO (e.g., NO x 
or Cl 2 or HC1) is added to the original fluid G. The concen- 
tration value c o (m0) of the selected molecule present in the 
original fluid G is unknown, and the configuration of the CNT 
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network is unknown. It is assumed that the response value 
difference AV increases approximately linearly with the con- 
centration difference Ac(mO) of the selected molecule so that) 


AU 1 (mO)=v 0 +v 1 -(c 0 (mO)+A 1 c(mO)), (Bl) 5 

AV 2 (mO)=v 0 +v l -(c 0 (mO)+A 2 c(mO)). (B2) 


A molecule mO should be chosen for which IA 2 c(mO)-A 1 c 
(mO)l is at least equal to a selected positive threshold. The 10 
quantities v 0 , v l and c 0 (mO) are then related by the equations 


vi = (AV2(*wO) - AV / i(mO))/{A2c(mO) - A ic(mO)}, (B3) 

v 0 + viCo(mO) = AV^mO) - viAic(mO) (B4) 

= AV 2 {mO) - viA 2 c(mO) 


and v 0 and v 1 are determined, in part, by the CNT network 20 
configuration (assumed fixed and reusable) that is present. 
Where, as is likely, v o i^0, one infers that the initial concen- 
tration value c o (m0) for the molecule mO is limited by 


c 0 (mO)={AU 1 (mO)-vO}/v 1 -A 1 c(mO)^{AF 1 (mO)/v 1 }- 25 

A^mO) (B5-1) 

or 

c 0 (mO)={A V 2 (mO)-vO }/v 1 -A 2 c(mO) ^{A V 2 (mO)/v 1 }- 30 

A 2 c(mO) (B5-2) 


Equations (B5-1) and (B5-2) provide an upper bound for the 
quantity c o (m0). Where it is known that the coefficient v 0 is 
non-positive, Eqs. (B5-1) and (B5-2) can be inverted to pro- 35 
vide lower bounds for the concentration: 

c 0 (mO)={A V l (mO)-vO }/v x - A { c(m0) = {A V { ( m0)/v l }- 

Ajc(m0) (B6-1) 


c 0 (mO)={A r 2 (m0)-v0 }lv i -A 2 c(m§) (A V 2 {m0)lv i }- 

A 2 c(w0) (B6-2) 

This Approach does not Provide a Direct Estimate for the 45 
Quantity C o (M0), Only an indication of whether the molecule 
mO is or is not present in a concentration of no more than the 
right hand quantity in Eqs. (B5-1) or (B5-2). However, this 
approach does not require determination and use of the 
response coefficients a z / that are required for the putatively 50 
more accurate method set forth in Appendix A. The method of 
Appendix B can be used to estimate upper (or lower) bounds 
for concentration values c of one, two or more selected mol- 
ecules. 

FIG. 7 is a flow chart of a procedure for estimating an upper 55 
bound for a concentration value c o (m0) for a selected mol- 
ecule in the fluid G. It is assumed that the response value 
difference AV varies approximately linearly with the concen- 
tration difference Ac(mO) of the selected molecule. In step 71, 
first and second (distinct) known increments, A^mO) and 60 
A 2 c(m0), of a selected molecule mO are added to a fluid G to 
provide first and second augmented fluids, G1 and G2. In step 
72, response value differences, 

A V { (mO))=V(Gl ;meas)- F(G;meas), (B7) 65 

(B8) 
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are measured or otherwise provided. In step 73, the coeffi- 
cient vl in an approximation for response value differences 


A V x (mO)=v 0 +v 1 ■{c o (m0)+A l c(m0) } , 

(B9) 

AV 2 (mO)=v 0 +v l -{c 0 (mO)+A 2 c(mO)}, 

(BIO) 

is determined according to 


v l =(AV 2 (m0)-AV l (m0))/{A 2 c(m0)-A 1 c(m0)}. 

(Bl 1) 


In step 74, the system queries whether the coefficient v 0 is 
likely non-negative. If the answer to the query in step 74 is 
“yes,” the system estimates an upper bound for the initial 
concentration value c o (m0), in step 75: 

c o (m0) =i{A V i (mO)/v 1 -A 1 c(mO)=AV 2 (wO)/v 1 -A 2 c(mO) 

(v o =0). (B12) 


If the answer to the query in step 74 is “no,” the system 
estimates a lower bound for the initial concentration value 
c o (m0), in step 76: 

c 0 (mO)^{AU 1 (wO)/v 1 -A 1 c(mO)=AV 2 (wO)/v 1 -A 2 c(mO) 

( v o=0). (B13) 

Appendix C. Effect of Varying Environmental Parameter. 

FIG. 18 is a flow chart of a procedure for estimating con- 
centration of a candidate fluid component, suspected to be 
present in a target fluid, where a response parameter V is 
associated with a controllably varied environmental param- 
eter (“EP”), and measurements of the response V(G;meas) are 
available for a sequence of reference values of EP, {EP^}^. 
The environmental parameter(s) EP may be one or more of 
the following: relative humidity RH, temperature T, gauge 
pressure GP, ambient chemicals present AC, and elapsed time 
At since the test sample was prepared. 

In step 191, measurements of the response V(EP^;meas) 
(p=0, 1 , . . . , P;Pi^ 1 ) are provided, where p=0 may correspond 
to EP equal to a reference value, such as RH=0 or 0.1, for 
example. In step 192, a sequence of reference values V(EP^; 
ref;h) (h=l , . . . , H:Hi^ 1) is provided, for each of H candidate 
fluid components and for each of P reference values of the 
environmental parameter EP. In step 193, the system deter- 
mines an error value 


2*W = {£ 

[P= 1 


V (EPp ; me as) - V (EP p ; ref-, h ) | 


<h=U ■■■ H), 


(Cl) 


where w p is a non-negative weight associated with the envi- 
ronmental parameter value EP^ and q is a selected positive 
number (e.g., q=l or 2). If desired, the approximation F e ((EP/ 
EP 0 )™;a) in Eq. (6) may be substituted for the reference quan- 
tity V(EP / ,;ref;h). 

In steps 194 and 195, a minimum error value (or values) 

e(min)=min{e(/z=l), . . . , e(h=H)} (C2) 

is/are determined (e.g., €(min)=€(h=h0)), and €(min) is com- 
pared with a selected threshold value e(min;thr). If 

e (min) =e (/z =/z 0) e (min ;thr) , (C3) 

the system interprets these conditions as indicating that the 
candidate fluid component (s) corresponding to h=hO is/are 
likely present in the target fluid, in step 206. If the condition 
(C3) is not satisfied, the system interprets this result as indi- 


A U 2 (wO))=U(G2;meas)- U(G;meas), 
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eating that none of the fluid component(s) hO is likely to be 
present in the target fluid, in step 197. 

What is claimed is: 

1 . A system for monitoring and evaluating status of an array 
of one or more underground mineral -producing wellsites, the 5 
system comprising: 

a first chip comprising: (i) a first nanosensor to estimate 
H 2 0 content present in the wellsite; and (ii) a second 
nano sensor to estimate presence or absence of each of at 
least first and second different target gas molecules in to 
the wellsite; 

a second chip containing (iii) a temperature microsensor to 
estimate local temperature in an underground fluid min- 
eral-producing wellsite; and at least one of: (iv) a pres- 
sure microsensor to estimate local fluid pressure adja- 15 
cent to the wellsite; (v) a distance microsensor to provide 
an estimate of distance to each of at least two spaced 
apart surfaces bounding the fluid mineral in the wellsite; 

(vi) a fluid flow microsensor to estimate at least one of 
flow direction and flow velocity of a fluid adj acent to the 20 
wellsite; and 

a wireless data transmission mechanism to receive and 
process data from the first and second nanosensors and 
from the at least one microsensor, and to transmit the 
processed data to a receiver spaced apart from the data 25 
transmission mechanism. 

2. The system of claim 1, wherein said second nanosensor 
is configured to estimate presence or absence of said at least 
first and second target gas molecules drawn from the group of 
carbon-containing molecules consisting of C m2 H w2 (m2i^l; 30 
n2=2m2, 2m2+l or 2m2+2), and C ml H wl O (ml 1 ; 
nl^2ml+2). 

3. The system of claim 1, wherein said first nanosensor 
comprises: 

an array of carbon nanosensor structures (CNSs), wherein 35 
at least one CNS in the array is loaded with a selected 
loading substance between first and second ends of the at 
least one CNS, and the first and second ends are con- 
nected: to at least one of a voltage source and a current 
source, and to an electrical parameter measurement 40 
mechanism that measures a change in at least one of 
electrical current, voltage difference, electrical resis- 
tance and electrical conductance between the first and 
second ends of the at least one CNS; and 

a computer programmed (1) to receive and compare the 45 
measured electrical parameter with a reference electrical 
parameter value and to estimate at least one of: (i) said 
H 2 0 content of said fluid adjacent to the at least one 
loaded CNS and (ii) presence or absence of each of said 
at least two target molecules in said wellsite; (2) to 50 
receive first and second incremental values, AV : and 
AV 2 , associated with a change in said measured electri- 
cal parameter associated with H 2 0 and with a represen- 
tative hydrocarbon, respectively; (3) to form a ratio AV X / 
AV 2 ; and (4) to interpret this ratio as approximately 55 
proportional to a/(l-a), where a (O^a^l) is a fraction 
of H 2 0 present in a mixture of fluids in a head space for 
said wellsite. 

4. The system of claim 3, wherein said at least one CNS is 
loaded with a coating on at least one surface of said at least 60 
one CNS. 

5. The system of claim 3, further comprising loading 
wherein said at least one CNS is loaded with a dopant distrib- 
uted within said at least one CNS. 

6. The system of claim 3, wherein said loading substance 65 
includes one or more of the substances hydroxypropyl cellu- 
lose and chlorosulfonated polyethylene. 
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7. The system of claim 3, wherein said loading substance 
includes one or more of the substances Pd, Pt, Ru, Rh, Ir, Os, 
Au, Ag and Hg. 

8. The system of claim 1, wherein at least one of said first 
target gas molecule and said second target gas molecule 
includes at least one of a hydrocarbon and an oxide of carbon. 

9. The system of claim 1, wherein said fluid flow micro sen- 
sor comprises first, second and third micro sensors, spaced 
apart and non-collinear, to sense at least one of (i) correlation 
of a flow of first and second fluid flow components and (ii) 
correlation of a flow of first and third fluid flow components. 

10. The system of claim 1, wherein said flow micro sensor 
is configured to measure at least one of said fluid flow direc- 
tion and said fluid flow velocity in response to injection of a 
pressurized fluid at a location spaced apart from said wellsite. 

11. The system of claim 1, wherein said second chip com- 
prises: 

an array of carbon nanosensor structures (CNSs), wherein 
at least one CNS in the array is loaded with a selected 
loading substance at one or more locations between first 
and second ends of the at least one CNS, and the first and 
second ends are connected to at least one of a voltage 
source and a current source, and are connected to an 
electrical parameter measurement mechanism that mea- 
sures a change in at least one of electrical current, volt- 
age difference, electrical resistance and electrical con- 
ductance between the first and second ends of the at least 
one CNS for each of said at least first and second target 
gas molecules; and 

a computer that is programmed to receive and compare the 
measured electrical parameter with a reference electrical 
parameter value and to estimate at least one of: relative 
humidity present at the at least one loaded CNS, tem- 
perature present at the at least one loaded CNS; fluid 
pressure present at the at least one loaded CNS. 

12. A method for analyzing status of an underground min- 
eral-producing wellsites, the method comprising: 

receiving first data at a first chip comprising: (i) a first 
nanosensor configured to estimate water content present 
in the wellsite; and (ii) a second nanosensor configured 
to estimate presence or absence of each of at least first 
and second different target gas molecules in the wellsite; 

receiving second data at a second chip containing at least 
one of: (iii) a temperature microsensor to estimate local 
temperature in an underground fluid mineral-producing 
wellsite; (iv) a pressure microsensor to estimate local 
fluid pressure adjacent to the wellsite; (v) a distance 
microsensor to provide an estimate of distance to each of 
at least two spaced apart surfaces bounding the fluid 
mineral; (vi) a fluid flow microsensor to estimate at least 
one of flow direction and flow velocity of a fluid adjacent 
to the wellsite; and 

processing the first data and the second data, and transmit- 
ting the processed data to a receiver spaced apart from 
the data transmission mechanism. 

13. The method of claim 12, further comprising drawing 
said target molecules from the group of molecules consisting 
of C m2 U„ 2 (m2i=l; n2=2m2, 2m2+l, 2m2+2), H 2 0 and 
C ml H wl O (mli=l; nl^2ml+2). 

14. The method of claim 12, further comprising providing 
said first nanosensor with an array of carbon nanosensor 
structures (CNSs), wherein at least one CNS in the array is 
loaded with a selected loading substance between first and 
second ends of the at least one CNS, and the first and second 
ends are connected: to at least one of a voltage source and a 
current source, and to an electrical parameter measurement 
mechanism that measures a change in at least one of electrical 
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current, voltage difference, electrical resistance and electrical 
conductance between the first and second ends of the at least 
one CNS; and 

providing a computer programmed to receive and compare 
the measured electrical parameter with a reference elec- 
trical parameter value and to estimate at least one of: (i) 
said water content of said fluid adjacent to the at least 
one loaded CNS and (ii) presence or absence of each of 
said at least first and second target gas molecules in said 
wellsite. 

15 . The method of claim 14 , wherein said computer is 
further programmed: 

to receive first and second incremental values, AV X and 
AV 2 , associated with a change in said measured electri- 
cal parameter associated with H 2 0 and with a represen- 
tative hydrocarbon, respectively; 

to form a ratio AV 1 /AV 2 ; and 

to interpret this ratio as approximately proportional to 
a/(l -a), where a (O^ia^i 1) is a fraction of H 2 0 present 
in a mixture of fluids in a head space for said wellsite. 

16 . The method of claim 14 , further comprising loading 
said at least one CNS by providing a coating of said loading 
substance on at least one surface of said at least one CNS. 

17 . The method of claim 14 , further comprising loading 
said at least one CNS with a dopant distributed within said at 
least one CNS. 

18 . The method of claim 14 , further comprising choosing 
said loading substance to include one or more of the sub- 
stances hydroxypropyl cellulose and chlorosulfonated poly- 
ethylene. 

19 . The method of claim 14 , further comprising choosing 
said loading substance to include one or more of the sub- 
stances Pd, Pt, Ru, Rh, Ir, Os, Au, Ag and Hg. 

20 . The method of claim 12 , wherein at least one of said 
first gas molecule and said second gas molecule is chosen 
from a group consisting of a hydrocarbon and an oxide of 
carbon. 

21 . The method of claim 12 , further comprising providing, 
as said flow microsensor, a mechanism comprising first, sec- 
ond and third microsensors, spaced apart and non-collinear, 
to sense at least one of (i) correlation of a flow of first and 
second fluid flow components, and (ii) correlation of a flow of 
second and third fluid flow components. 

22 . The method of claim 12 , further comprising using said 
flow micro sensor to measure at least one of said fluid flow 
direction and said fluid flow velocity in response to injection 
of a pressurized fluid at a location spaced apart from said 
wellsite. 

23 . The method of claim 12 , further comprising providing, 
as said second chip, a mechanism comprising: 

an array of carbon nanosensor structures (CNSs), wherein 
at least one CNS in the array is loaded with a selected 
loading substance at one or more locations between first 
and second ends of the at least one CNS, and the first and 
second ends are connected to at least one of a voltage 
source and a current source, and are connected to an 
electrical parameter measurement mechanism that mea- 
sures a change in at least one of electrical current, volt- 
age difference, electrical resistance and electrical con- 
ductance between the first and second ends of the at least 
one CNS for each of said at least two target molecules; 
and 

a computer that is programmed to receive and compare the 
measured electrical parameter with a reference electrical 
parameter value and to estimate at least one of: relative 
humidity present at the at least one loaded CNS, tern- 
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perature present at the at least one loaded CNS; fluid 
pressure present at the at least one loaded CNS. 

24 . The method of claim 12 , wherein said step of providing 
said second nanosensor comprises: 

5 providing first and second interdigitated electrodes, which 
are connected to at least one of (i) a voltage source 
having a controllable voltage difference and (ii) a con- 
trollable current source; 

providing a zeroth order electrically conducting network, 
to numbered 0, of uncoated single wall carbon nanotubes 
(“SWCNTs”) forming a path P0 between the first and 
second electrodes; 

providing N electrically conducting networks, numbered 
n=l, . . . , N(Ni^l) of single wall carbon nanotubes 
1 5 (“SWCNTs”) forming paths Pn (n=l, . . . , N) between 
the first and second electrodes, where at least one 
SWCNT in each of the networks no. n (n=l, . . . , N) is 
loaded with molecules of a selected chemical no. n; 

providing in the chamber a gas G that may have molecules 
20 of a gas constituent no. m present in a concentration C m 
(m=l, . . . , M; Mi^l), where at least one of the concen- 
tration values C m is unknown, allowing at least one mol- 
ecule of the gas G to become absorbed on at least one of 
the SWCNTs in each of the paths P0 in the zeroth order 
25 network and Pn in the network number n (n=l, . . . , N); 

providing an electrical parameter value (“response value”) 
V 0 (G;meas), associated with the path P0, of a change in 
at least one of (i) electrical current, (ii) electrical con- 
ductance, (iii) voltage difference and (iv) electrical resis- 
tance associated with the path P0, when the gas G is 
present, and providing a corresponding response value 
V„(G;meas), associated with the path Pn (n=l, . . . , N) 
when the gas G is present; 

adding a known concentration increment A x C m0 of a gas 
constituent no . mO , selected from among the M constitu- 
ents, to the gas G to provide an augmented gas G', 
allowing at least one molecule of the gas G' to become 
absorbed on at least one of the SWCNTs in each of the 
path P0 and Pn (n=l, . . . , N), and providing correspond- 
40 ing response values V 0 (G';meas) and V w (G';meas) for 
the gas G'; and 

comparing differences, V w (G;meas)-V 0 (G;meas) and 
V w (G;meas)-V 0 (G , ;meas) (n=l, . . . , N) of the response 
values and estimating at least one of (i) concentration 
45 and (ii) concentration upper bound and (iii) concentra- 
tion lower bound for at least one of the first gas molecule 
and the second gas molecule present in at least one of the 
gas G and the gas G\ 

25 . The method of claim 24 , wherein said process of esti- 
50 mating at least one of (i) said concentration C and (ii) said 

concentration upper bound and (iii) said concentration lower 
bound comprises: 

providing an error function c(C l5 . . . , C M ), defined by 

55 2e(C x , . . . , C M )=2„vv„d^,(G;meas)-F 0 (G;meas)- 

2 m a n ra C ra } 2 +2„w , „-{ r„(G';meas)- r 0 (G';meas)- 
a n^(C m + ACJ} 2 , 

where a n m is a response value coefficient relating change in 
said response value V„(G;meas) to change in said concentra- 
60 tion C m of said constituent no. m, AC m is a known concentra- 
tion increment of said gas constituent no. m, added to said gas 
Gto provide said gas G’, w„, and w’„ are selected non-negative 
weight values, and at least one increment value AC m is posi- 
tive; 

65 providing non-negative values, C m (min) for said concen- 
trations C m that minimize the value of the error function 
e(C L , C M ); and 
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associating at least one of the values C m (min) with an 
estimate of said concentration of said constituent no. m 
in at least one of said gas G and said gas G\ 

26. The method of claim 25, further comprising choosing 
said weight numbers w w and w’ w to be substantially equal for 
at least one of said integers n. 

27. The method of claim 24, wherein said process of esti- 
mating at least one of (i) said concentration C and (ii) said 
concentration upper bound and (iii) said concentration lower 
bound comprises: 

adding a second known concentration increment A 2 C m0 of 
said selected gas constituent, selected from among the M 
constituents, to said gas G to provide a second aug- 
mented gas G", allowing at least one molecule of the gas 
G" to become absorbed on at least one of the SWCNTs 
in each of the paths PO and Pn (n=l, . . . , N), and 
providing corresponding response values V 0 (G M ;meas) 
and V w (G";meas) for the gas G"; 

providing response value differences, A V w (G ; G)=V W (G ; 
meas)-V„(G;meas) and AV w (G M, ;G)=Y K (G M ;meas)-V w 
(G;meas); 

determining a response coefficient vl as 

vl ={AV(G";G)-AV(G';G)}/{A 2 C m 0 -A 1 C m0 }; and 

estimating said concentration C m0 of said selected gas con- 
stituent no. mO in said gas G according to at least one 
relation 

C wl0 ={Ar(GVG)-vO}/vl-A 1 C w , 0 } 

and 

C m 0 ={AV(G";G)-rt}/vl-A 2 C m0 }, 

where vO is an estimated reference value of a response value 
difference, AV„(G';G) or AV W (G";G). 

28. The method of claim 27, further comprising estimating 
an upper bound for said concentration CmO according to at 
least one relation 

C m 0 ^{AV(G';G)ivl}-A L C m0 

and 

C m o={AV(G";G)/vl}-A 2 C m0 , 

where said estimated reference value vO is non-negative. 

29. The method of claim 27, further comprising estimating 
an upper bound for said concentration CmO according to at 
least one relation 

GwO= { A V(G f ; G)/vl }A l C m } 

and 

GhO={A V(G”;G)lv\ }-A 2 C m0 , 

where said estimated reference value vO is non-positive. 

30. A system for monitoring and evaluating status of an 
array of one or more underground mineral-producing 
wellsites, the system comprising: 

a first chip comprising: (i) a first nanosensor to estimate 
H 2 0 content present in the wellsite; and (ii) a second 
nano sensor to estimate presence or absence of each of at 
least first and second different target gas molecules in 
the wellsite; 
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a second chip containing (iii) a temperature micro sensor to 
estimate local temperature in an underground fluid min- 
eral-producing wellsite; and at least one of: (iv) a pres- 
sure microsensor to estimate local fluid pressure adja- 
5 cent to the wellsite; (v) a distance microsensor to provide 

an estimate of distance to each of at least two spaced 
apart surfaces bounding the fluid mineral in the wellsite; 
(vi) a fluid flow microsensor to estimate at least one of 
flow direction and flow velocity of a fluid adjacent to the 
to wellsite; 

where the fluid flow microsensor comprises first, second 
and third microsensors, spaced apart and non-collinear, 
to sense at least one of (i) correlation of a flow of first and 
second fluid flow components and (ii) correlation of a 
15 flow of first and third fluid flow components; and 

a wireless data transmission mechanism to receive and 
process data from the first and second nanosensors and 
from the at least one micro sensor, and to transmit the 
processed data to a receiver spaced apart from the data 
20 transmission mechanism. 

31. A system for monitoring and evaluating status of an 
array of one or more underground mineral-producing 
wellsites, the system comprising: 

a first chip comprising: (i) a first nano sensor to estimate 
25 H 2 0 content present in the wellsite; and (ii) a second 

nano sensor to estimate presence or absence of each of at 
least first and second different target gas molecules in 
the wellsite; 

a second chip containing (iii) a temperature micro sensor to 
30 estimate local temperature in an underground fluid min- 
eral-producing wellsite; and at least one of: (iv) a pres- 

sure microsensor to estimate local fluid pressure adja- 
cent to the wellsite; (v) a distance microsensor to provide 
an estimate of distance to each of at least two spaced 
35 apart surfaces bounding the fluid mineral in the wellsite; 
(vi) a fluid flow microsensor to estimate at least one of 
flow direction and flow velocity of a fluid adjacent to the 
wellsite; 

where the second chip further comprises an array of carbon 
40 nanosensor structures (CNSs), wherein at least one CNS 
in the array is loaded with a selected loading substance at 
one or more locations between first and second ends of 
the at least one CNS, and the first and second ends are 
connected to at least one of a voltage source and a 
45 current source, and are connected to an electrical param- 

eter measurement mechanism that measures a change in 
at least one of electrical current, voltage difference, elec- 
trical resistance and electrical conductance between the 
first and second ends of the at least one CNS for each of 
50 said at least first and second target gas molecules; 

a computer that is programmed to receive and compare the 
measured electrical parameter with a reference electrical 
parameter value and to estimate at least one of: relative 
humidity present at the at least one loaded CNS, tem- 
55 perature present at the at least one loaded CNS; fluid 
pressure present at the at least one loaded CNS; and 
a wireless data transmission mechanism to receive and 
process data from the first and second nanosensors and 
from the at least one micro sensor, and to transmit the 
60 processed data to a receiver spaced apart from the data 
transmission mechanism. 



